Induced magnetoelectric response in Pnma perovkites 
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We use symmetry analysis to show that the G, C and yl-type antiferromagnetic Pnma perovskites 
can exhibit magnetoelectric (ME) responses when a ferroelectric instability is induced with epitaxial 
strain. Using first-principles calculations we compute the values of the allowed ME response in 
strained CaMnOa as a model system. Our results show that large linear and non-linear ME responses 
are present and can diverge when close to the ferroelectric phase transition. By decomposing the 
electronic and ionic contributions, we explore the detailed mechanism of the ME response. 



Interest in magnetoelectric (ME) materials has in- 
creased over the last few years because of their cross 
coupling between the electric polarization and magne- 
tization and their consequent potential for technological 
applications [1]. However, the search for good MEs is 
facing difficulties: compounds with the required sym- 
metry (breaking of both the time and space inversion) 
are uncommon, and when these requirements are met, 
it is often at low temperatures. In addition, although 
the magnitude of the response is in principle bounded by 
the product of the dielectric and magnetic permeabilities 
(a < a/c/x), in practice it tends to be much smaller than 
this value. One promising direction in the search for im- 
proved magnetoelectrics is the exploration of multiferroic 
materials, since the presence of multiple ferroic orders 
often presents the desired coupling properties for large 
ME responses [1-3]. Another route is the engineering of 
artificial heterostructures with specific chemistries and 
symmetries [4-6]. Here we demonstrate from symmetry 
considerations that the Pnma G, C or ^4-type antiferro- 
magnetic perovskites, which are not multiferroic and do 
not allow a ME response in their bulk form, can become 
ME when a polar distortion is induced using thin film 
heteroepitaxial strain. Then, using first-principles calcu- 
lations for a model example - Pnma CaMnOa - we show 
that particularly large ME responses can be achieved in 
the vicinity of the ferroelectric phase transition. 

Technical details — We performed all calculations 
within density functional theory as implemented in the 
VASP code [7, 8]. Since the purpose of the present study 
is to provide a model example, we restricted ourselves 
to the Local Density Approximation (LDA) functional 
and intentionally avoided studying the U and J depen- 
dence of the LDA+C7 functional [9]. We oriented the 
Pnma unit cell with the longest axis along the b direc- 
tion, and applied a cubic epitaxial strain on the a and 
c directions by imposing a — c and relaxing the b di- 
rection. To compare with previous studies on strained 
CaMnOa [10], we performed all calculations at the calcu- 
lated generalised gradient approximation (GGA) PBEsol 
volumes for each strain (see Supplemental Information). 
The non-collinear properties and the ME response were 
well converged at a plane wave cutoff of 550 eV and a 



4x2x4 k-point grid. To compute the ME responses, we 
applied a Zeeman magnetic field to the spins with the 
spin-orbit coupling included and extracted the polariza- 
tion response as described in Ref. 11. 

Group theory analysis — The Pnma structure is 
obtained by condensation of three antiferrodistortive 
(AFD) instabilities coming from the M and R zone 
boundary points of the high symmetry ^4BX 3 cubic cell 
[12]; The primitive unit cell contains four formula units. 
When considering only one type of magnetic atom at the 
Wykoff site 4a in the Pnma cell, four collinear magnetic 
orders are likely: ferromagnetic (F) and three antiferro- 
magnetic (AFM) of G, C and A types. For each collinear 
magnetic order we consider spin directions along x, y 
or z, corresponding to the a, b and c axes of the crys- 
tal. This allows us to define 3x4 magnetic order pa- 
rameters. In Tab. I we report the symmetry characters 
of each of these magnetic order parameters according to 
their transformation under the symmetries of the Pnma 
space group. Interestingly, the magnetic order param- 
eters group in threes with the same character. In an 
expansion to second order of the energy with respect to 
these order parameters, the couplings between the three 
order parameters with the same character are allowed. 
According to Tab. I, we then conclude that for G, C, A 
or F magnetic order in the Pnma space group and for 
any easy axis, spin cantings in the two other directions 
are allowed. We note that if we consider only the AFM 
magnetic orders, then three characters out of four allow 
F spin canting i.e. for weak ferromagnetism (FM). 
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TABLE I. Symmetry character of the magnetic order param- 
eters of the perovskite 20 atom unit cell with the Pnma and 
Pmc2\ space groups following the standard settings as given 
in the Bilbao Crystallographic Server. When it is allowed, 
we report the components of the linear ME tensor that are 
non-zero. 

As shown in Tab. I, if we consider the condensation 
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of any of the magnetic order parameters in the Prima 
structure, none of the resulting magnetic space groups 
permits for ME response. This is because, although the 
time-reversal symmetry is broken, the inversion center 
is still preserved. One possibility to break the inversion 
symmetry is with polar displacements. Recently, it was 
shown from first-principles that it is possible to induce 
a ferroelectric (FE) instability in Prima CaMnC>3 with 
epitaxial strain [10]. When the polarization develops, 
the system condenses into the Pmc2\ space group which 
has no inversion center. In Tab. I we report the symme- 
try characters of the magnetic order parameters in the 
Pmc2\ space group. Interestingly, the magnetic orders 
stay in the same groups of three as in the Prima struc- 
ture. Three out of four characters in the new Pmc2\ 
space group allow linear ME coupling. This means that 
in G, G and A-type AFM Prima perovskites with one 
type of magnetic cation, it is possible to induce a ME re- 
sponse if a FE instability develops with epitaxial strain. 

Ground state properties — Most previous studies on 
CaMnC>3 did not explore the possibility of non-collincar 
spin canting; to our knowledge there is only one old ex- 
perimental report of weak FM [13]. Therefore we first 
performed calculations for bulk CaMnC>3 including the 
spin-orbit coupling and explicitly checking all possible 
collinear and non-collinear magnetic ground states. We 
found that the lowest energy magnetic ordering is G-type 
AFM with the easy axis along the z direction. This in- 
deed allows a canting of the spins of the F-type along 
the y direction and of the A-type AFM along the x direc- 
tion (G z F y A x ground state), consistent with our group 
theory analysis reported in Tab. I and the experimen- 
tal measurements [13]. The calculated canting angles are 
2.6° along the x direction and 1.0° along z. resulting in 
a weak FM of 0.04 fj,B per Mn atom. This calculated 
weak FM is larger than the experimental report (0.004 
Ms [13]); the discrepancy could be due to experimental 
uncertainty or our use of the LDA approximation. 

Having verified the accuracy of our calculated ground 
state properties, we then checked how these properties 
are affected by the epitaxial strain and by the devel- 
opment of the FE polarization. To determine the low- 
est energy phase, we performed full atomic and spin re- 
laxations at different epitaxial strains and for all of the 
magnetic order parameters reported in Tab. I. We sum- 
marize the results in the phase diagram of Fig. 1. As 
reported previously [10], we oberve that a FE instability 
appears beyond critical tensile epitaxial strain rj FE . r\ FE 
was predicted previously within the GGA Wu-Cohen ap- 
proximation to be 2.0%, the LDA value obtained here is 
larger (3.2%). This shows that the magnitude of strain 
required to induce ferroelectricty is strongly dependent 
on the approximations used (exchange-correlation func- 
tional, plane wave scheme, etc.) in a calculation, and ex- 
act quantitative predictions should be made with caution. 
The electric polarization develops along the in-plane c 



direction and increases with 77, reaching large values at 
large strains (27 /zC.cm -2 at 77=4.5%). Interestingly, we 
found that the amplitudes of the AFD rotations are only 
slightly modified (few %) by the epitaxial strain and per- 
sist even when the polarization develops. 

Looking at the magnetic properties we found that the 
strain induces a spin flop of the easy axis from the z di- 
rection to the y direction at a critical strain rff = 2.6% 
(Fig. 1) smaller than r/ FE . This transition keeps the pri- 
mary G-type AFM magnetic order but, since the easy 
axis is modified, changes the type of spin canting allowed 
from G z F y A x to G y C x F z . With G y order, canting of the 
G-type along the x direction and a weak FM along the z 
directions are allowed (Tab. I). When the FE transition 
takes place, we find no change of the magnetic ground 
state, with the system staying in the G V C X F Z magnetic 
state. We observe however that as the polarization in- 
creases, the total magnetic moment and spin cantings 
decrease by respectively a few % and a factor of three 
(between 77 = 3.2% and 4.5%). We confirm that this is 
an effect of the FE polarization and not of the strain by 
constraining the polarization to be zero by symmetry - 
in this case we do not observe such a reduction of the 
spin cantings. At large tensile epitaxial strain, we ob- 
serve a magnetic phase transition from the G-type AFM 
to the A-type AFM with again a change of the easy axis 
from the y (G y ) to the z direction (A z ) at r) GA = 4.6%. 
Here again, consistent with Tab. I, with the A z order we 
observe spin cantings of the G-type along the y direction 
and of the G-type along the x direction (A z C y G x ). The 
weak FM is then lost during this phase transition. 
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FIG. 1. Multiferroic phase diagram of CaMnOa under tensile 
epitaxial strain. The top part shows the magnetic orders and 
the lower part shows the crystallographic orders (PE = para- 
electric). rj sf = 2.6%, n FE = 3.2% and rf A = 4.6%. Between 
77 — and 77 s ' the magnetic point group is m'mm' , between 
Tff and n FE it is m'm'm, between rj FE and n GA it is m'm'2 
and beyond n GA it is mm2. 

ME response — From the phase diagram of Fig. 1 and 
according to the group theory analysis of Tab. I we can 
see that in the region between r\ FE and rj GA the sys- 
tem allows for a diagonal linear ME response [G y C x F z 
magnetic state in the Pmc2\ crystal space group). How- 
ever, such analysis does not allow us to predict the am- 
plitude of the response. To determine the amplitude of 
the ME response, we performed calculations under a Zee- 
man magnetic field along the x, y and z directions and 
calculated the induced electric polarization [11]. In Fig. 2 
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we report the induced change in polarization versus mag- 
netic field at three different epitaxial strains between rj FE 
and i] GA . Fig. 2. a shows the variation of the polarization 
(AP) along the x direction (AP X = P x since P x = at 
B = 0) when a magnetic field is applied in the same direc- 
tion, giving the a xx component of the ME response. We 
find a linear evolution of the polarization for B x between 
-25 and 25 T and ME response values: a xx — —16 ps.m -1 
at i] = 3.3% and —12 ps.m -1 at -q = 4.5%. These values 
are large compared with the prototypical ME compound 
Cr2C>3 where the calculated ME response is 1.45 ps.m -1 
[11]. In addition to the polarization induced along the x 
direction, we also observe a change in the FE polarization 
along the z direction even if the field is applied along the 
x direction. We report this response in Fig. 2. a where 
we see that it is highly non- linear. This is in agreement 
with group theory wich shows that linear response a xz 
is not permitted, but that the next order non- linear re- 
sponse (3 XXZ [14] is allowed. This non- linear response is 
small with respect to the linear response and is extremely 
sensitive to the epitaxial strain (j3 xxz = -7xl0 -19 s/A at 
3.3%). In Fig. 2.b we report the variation of polarization 
along the z direction when a magnetic field is applied in 
the same direction (a zz component). Here again the po- 
larization response is strongly sensitive to the epitaxial 
strain and deviates from an ideal linear response, con- 
sistent with symmetry analysis which yields a non-zero 
a zz and fi zzz . The a zz value is particularly large near to 
the FE transition (a zz — —85 ps.m -1 at r\ = 3.3%, close 
to rj FE ) and decays rapidly away from r/ FE (a zz = —19 
ps.m -1 at r] = 4.5%) and (3 ZZZ is sizeable (—4.2 x 10 -16 
s/A at 3.3%). Finally, when applying a magnetic field 
along the y direction, we do not observe a polarization 
along the y direction. This shows that even when a com- 
ponent is symmetry allowed, it can be extremely small in 
amplitude, here lower than the precision of our calcula- 
tions. With the field applied along y, we observe however 
a tiny polarization response along the z direction which 
corresponds to the allowed non-linear ME tensor compo- 
nent (3yy Z . 

Electronic versus ionic contribution — Recently it has 
been shown by explicit calculation of the Zccman elec- 
tronic contribution to the ME response (a elec ) in Cr 2 03 
and LiNiPC>4, that a elec can be comparable in magnitude 
to the ionic contribution (a %ontc ) [11]. It has also been 
suggested that in FE perovskites such as BiFeC>3 , the ME 
response is dominated by the FE soft mode (ionic contri- 
bution) [2, 3, 15]. To check how large are the electronic 
and ionic contributions to the ME response in strained 
CaMnOa, we computed these two contributions as in 
Ref.[ll]. We summarize the results in Fig. 3. a and b, 
where the ionic and electronic contributions to a xx and 
a zz are plotted with respect to the epitaxial strain. We 
can conclude the following: First, the ionic contribution 
clearly dominates the total response. a?,° mc is extremely 
sensitive to the strain and diverges when approaching 
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FIG. 2. Polarization response of CaMnOa under magnetic 
field at three different epitaxial strains, (a) Variation of po- 
larization along the x direction (empty symbols, plain lines) 
and along the z direction (plain symbols, dashed lines) when 
the field is applied along x and (b) Variation of polarization 
along the z direction when the field is applied along z. 

the FE phase transition, consistent with the softening of 
a polar mode along the z direction [3]. a x ° x tc is how- 
ever much less sensitive to the strain and does not show 
any divergence when approaching the FE transition. We 
find that the response of a xx is mainly dominated by 
a relatively soft mode (110 cm -1 ) which keeps almost 
the same frequency for epitaxial strain from r] FE to r\ GA . 
Second, the electronic contribution (Fig. 3.b) shows the 
opposite trend to the ionic contribution in both a xx and 
a zz . While the ionic contribution has the tendency to 
decrease when the strain increases, the electronic con- 
tribution increases. We remark that, even though the 
electronic contribution is much smaller than the ionic 
contribution, its absolute value is large compared with 
the values reported for Cr 2 03 (a elec — 0.34 ps.m -1 ) and 
LiNiP0 4 (a elec = 1.1 ps.m -1 ) [11]. In strained CaMn0 3 
the electronic contribution alone can even be larger than 
the total ME response of Cr 2 3 (a tot = 1.45 ps.m -1 
[11])- 

The opposite trend between a lomc and a elec shows 
clearly that the electronic response is not driven by a 
phonon instability. We can make an analogy with the 
dielectric permittivity in crystals (e r ). The ionic contri- 
bution to e is directly affected by the softening of polar 
phonon modes and is responsible for the divergence of e r 
close to a FE phase transition. However the electronic 
contribution (or high frequency e£°) does not diverge at a 
FE phase transition, but rather diverges close to an elec- 
tronic instability such as a metal-insulator (M-I) phase 
transition. We expect similar behaviour for a elec , how- 
ever we are unable to test it in strained CaMnOa as it 
does not show a M-I transition. We can, nevertheless, 
comment on the dependence of a elec on the gap. The 
magnetic phase transition between G y C x F z and A z C y G x 
at i] GA — 4.6% is first order and ji, a elec, s and ef ec, s 
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FIG. 3. xx (squares) and zz (circles) components of the ME 
(plain symbols, in ps/m) and dielectric constant (empty sym- 
bols, SI units) versus epitaxial strain: (a) ionic contribution, 
(b) electronic contribution, (c) magnetic permeability (tri- 
angles, SI units) and ME upper bound ^/ejl (stars, in ns/m 

whith e — £o(£r° n + £r' ec )) versus epitaxial strain, (d) elec- 
tronic bang gap (eV). 

reported in Fig. 3.b do not show diverging behaviour. 
While the decrease of is directly related to the in- 
crease of the band gap (see Fig.3.d), a elec in fact increases 
with decreasing gap (See Fig. 3.b). The magnetic perme- 
ability /i (we have here \x mx = \i zz = ^ elec = ^) decreases 
in the same range of strain (Fig3.c). This shows that 
the link between a elec and e^°, /i and the band gap is 
not straightforward. While a formulation of the orbital 
contribution to a [16] and a generalized Lyddane-Sachs- 
Teller relationship between a, /i and e in ME [17] have 
been reported, an exact formulation of Zeeman a elec is 
still needed. In Fig.3.c we also report the upper bound 
of a which is equal to yfeji [18] . As we can see ^fe\x can 
be as large as 50 ns/m while a is of the order of 90 ps/m. 
a is then about 0.1 % of its upper bound as reported for 
Cr 2 3 [18]. 

From these observations we propose the following de- 
sign rules to obtain large ME responses: (i) one or more 
soft polar modes to increase a through e* on , (ii) a mag- 
netic instability to increase a through fi (iii) proximity 
to an electronic instability such as a M-I phase transi- 
tion to increase e 1 ^ . (i) is clearly evidenced in our results 
on strained CaMnOa. (ii) is not observed in spite of the 
magnetic phase transition at -q GA . This is because the 
magnetic phase transition at rj GA is of the first order 
which does not result to the divergence of ix. However, 
second order magnetic phase transitions should have the 
diverging effect on fi and a. (iii) will cause to diverge 
but we cannot yet conclude how it will affect a and fur- 



ther explorations are needed to clarify this point. If (i), 
(ii) and (iii) can be achieved simultancoulsy, a very large 
ME response could be obtained if the effects on a are 
of the same sign. While in this paper we have shown 
the divergence of the ME response with epitaxial strain, 
the effect will also occur at phase transitions induced by 
temperature, pressure, etc. 

We emphasise that our findings are not restricted to 
the case of CaMnC>3, but are valid for all G, C and ^4-type 
AFM Prima perovskites since the group theory analysis 
reported in Tab. I is valid for any A. B and X ions. Since 
the Pnma structure is the most common natural ground 
state of the ABX3 perovskites [19], this offers many pos- 
sibilities for creating new ME materials with epitaxial 
strain. Furthermore, mixing chemistries in superlattices 
vastly increases the possibility of generating phases with 
coexisting phonon, electronic and magnetic instabilities 
and giant ME responses [4-6]. 
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